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1 Introduction 



The recent development 0-0 in the theory of integrable many-body systems is mainly 
related with the discovery || of dynamical r-matrices, i.e. r-matrices depending on phase 
variables. One natural way to understand the origin of dynamical r-matrices is to con- 
sider the reduction procedure [^D[-[P^1[,H, [|. In this approach one starts with an initial 
phase space V supplied with a symplectic action of some symmetry group. Considering 
a relatively simple invariant Hamiltonian and factorizing the corresponding dynamics by 
the symmetry group one gets a smaller phase space V re d with a nontrivial dynamics. 
Then the L-operator coming in the Lax representation ^ = [M, L] appears as a specific 
coordinate on V re d while the dynamical r-matrix describes the Poisson (Dirac) bracket on 
the reduced space. 

At present the reduction procedure is elaborated for the majority of integrable many- 
body systems and the corresponding r-matrices are derived. One of the most interesting 
exception is the elliptic Ruijsenaars- Schneider model [HJ. Recently two different dynam- 
ical r-matrices for this model were found in and ||. Both of these r-matrices were 
obtained by a direct calculation and the question of their equivalence still remains open. 

In this letter we apply the Poisson reduction procedure to the affine Heisenberg Dou- 
ble (HD) IJTBJ and derive the elliptic Ruijsenaars-Schneider model with the dynamical 
r-matrix. The reason to use the affine HD becomes apparent due to its relation with inte- 



grable many-body systems of Calogero type. As was shown in [17], [H| the Calogero-Moser 
and the rational and trigonometric Ruijsenaars-Schneider hierarchies can be obtained by 
means of the reduction procedure from the cotangent bundle of an affine Lie group T*G(z) 
and from a finite dimensional Heisenberg Double. The affine Heisenberg Double may be 
regarded as a deformation of T*G(z) and therefore one can suggest that the affine HD is a 
natural candidate for the phase space standing behind the elliptic Ruijsenaars-Schneider 
system. 

The plan of the paper is as follows. In the second section we briefly describe the affine 
HD in terms of variables which are suitable for the reduction procedure. Then we fix 
the momentum map, corresponding to the natural action of the affine Poisson-Lie group 
on HD. The solution of the momentum map equation is shown to be equivalent to the 
L-operator of the elliptic Ruijsenaars-Schneider model. In the third section we study the 
Poisson structure of the reduced phase space and prove that it coincides with the one 
of the elliptic Ruijsenaars-Schneider model. The dynamical r-matrix naturally appears 
in our consideration and is equivalent to the one obtained in 0. In the last section we 
show that the problem of solving the equations of motion is equivalent to the specific 
factorization problem. 

In our presentation we omit the detailed description of the HD and the proof of some 
statements. The complete discussion will be given in a forthcoming publication. 

2 Affine Heisenberg Double 

The general construction of a Poisson manifold known as the Heisenberg double was 



elaborated in [R]]. We shall discuss the HD for the affine GL(N). It is convenient to 
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describe the Poisson structure of the affine HD in the following form. Let A(x) and C(x) 
be formal Fourier series in a variable x with values in GL(N). The matrix elements of 
the harmonics of A(x) and C(x) can be regarded as generators of the algebra of functions 
on HD. The Poisson structure on HD looks as follows: 

-{A 1 (x),A 2 (y)} = -r T (x-y)A 1 {x)A 2 (y)-A 1 {x)A 2 {y)r ± {x-y) 



+ A 2 (y)r + (x -y- 2A)A 1 (x) + A 1 (x)r^(x -y + 2A)A 2 {y), 
{C7i(x),C 2 (j/)} = -r T (x-y)C 1 (x)C 2 (y)-C 1 (x)C 2 (y)r ± (x-y) 
+ C 2 (y)r + (x - y)Ci(x) + C x (x)r-{x - y)C 2 (y), 

-{A x (x),C 2 (y)} = -r_{x-y)A 1 {x)C 2 {y)-A 1 {x)C 2 {y)r„{x-y + 2A) 

7 

+ C 2 (y)r+(x - y)A 1 (x) + A x {x)r„{x -y + 2A)C 2 (y), 

-{d(x),A 2 (y)} = -r + {x-y)C l (x)A 2 (y)-C 1 (x)A 2 (y)r + {x-y-2A) 

7 

+ A 2 (y)r + (x -y- 2A)C 1 {x) + d(x)r_(z - y)A 2 (y), 

where 7 and A are complex numbers with Im A > 0. Here we use a standard tensor 
notation. The matrices r±(x) are defined by their Fourier series as 

r+ ( x ) = r+ + Pj2 e~ inx , r_(ar) = r_ - P ^ e inx , 

n>0 n>0 

where ^ 

r+ = - E u ® + ZJ E v ® 

j i<jr' 

r_ = —Pr + P and P is the permutation operator. It can be easily checked that 

r + (x) — r_(x) = 27rP5(a;) and Pr + (x)P = — r_(— x). 
In the region of convergence r±(x) coincide with the standard trigonometric r-matrix for 



the affine Lie algebra. The Poisson subalgebra generated by A(x) was introduced in |L9 

— * 

to describe the Poisson structure of GL(N) . 
Assuming the expansions 

A(x) = I + 7 J(x) + ■ ■ ■ , C(x) = g(x) + 2A = 7 fc, 

where A; is a (fixed) central charge, in the deformation limit 7 — > we recover the standard 
Poisson structure on the cotangent bundle T*GL(N) over the level k centrally extended 
current group GL(N). 

The action of the current group GL(N) on HD: 

A(x) -► T-\x - A)A(x)T(x + A), 
C{x) -> T-\x- A)C{x)T{x- A) 

is Poissonian. Thereby, we can consider the Poisson reduction of HD over the action of 
GL(N). 
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The momentum map taking value in GL(N) reads as follows: 

M{x) = A'\x - A)C(x - A)A(x - A)C~ 1 (x + A). 

It is easy to check that M(x) does generate the action of the current group. We fix the 
value of M(x) as: 

M(x) = e ih {l - 2m5 £ {x) 1 ~ e - K). (2.1) 

% 

Here h and e are arbitrary complex numbers, 

1 f - _ x 1 n=+oo i 

Se(x) = - (0(x + -)-9(x--) )= — £ -(e m f-e-^)e- 



2' ^ 2 V 2me n t- 00 n 

and is a constant matrix K = e ® e*, where e is the A^-dimensional vector with entries 
ei = l/ViV0 

Although eq.( |2.1|) can be solved for any value of e, the reduced phase space remains to 
be infinite dimensional after performing the factorization procedure . To extract a finite 
dimensional phase space let us carry out the following trick. By multiplying the both 
sides of fl2.1p on C(x + A), one gets 

i _ e ~ ix 

C(x + A) - e- ih A-\x - A)C{x - A)A{x - A) = 2m6 e {x)K : C{x + A). (2.2) 

The l.h.s. of this equation does not have any explicit dependence on e. As to the r.h.s., 
when e tends to zero, 5 e (x) tends to S(x) and the r.h.s. is well defined only if the 
function 1 ~ e . — C(x + A) is well defined at x = 0. Hence, this equation can be solved 
only for meromorphic functions C[x + A) with poles of the first order. In this case 
lim £ ^o 5 £ (x)±=f^C(x + A) = 8(x)Res x = C(x + A). 

So we define the constraint surface as being the solution of the equation 

C(x + A) - e- ih A~\x - A)C{x - A)A{x - A) = 2m5(x)KRes x=0 C(x + A) 

and in the following we shall explore solutions of this equation. 
We start with the following difference equation 

C{x + A) - e- ih D- 1 G{x - A)D = 2ni5{x)Y, (2.3) 

where D is a constant diagonal matrix and Y is an arbitrary constant matrix. Performing 
the Fourier expansion we get a solution of fl2.3|) in the form 

n=oo ginx 
ij n=-oo e e e 

where we use the notation Sy = h + qij, D = e tq , = q% — qj- 

1 It is worthwhile to note that in the deformation limit 7^0, - — > const, — — > const the constraint 



(2.1) reduces to the one used in pTiJ to get the trigonometric Ruijsenaars model. 
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It is useful to introduce the function of two complex variables 



w(x, s) 



n=co inx 

; y . 

pinA p — iSf,— inA 

n=-oo c c c 



It is clear that w(x, s) is a meromorphic function of s for any x : |Im x\ < Im A and has 
two obvious properties 

1. w(x, s + 2n) = w(x, s), 

2. w(x, s + 2A) = e tA - ix w{x, s). 

Moreover, as a function of s it has simple poles at 0, ±2ir, ±4tt, . . . and ±2A, ±4A, . . ., 
Res s= ow = 1. By these data w is uniquely defined as: 

, . a(s + x-A) _cW/ T _A t , fn a\ 

w(x,s) = -j . , . ; e - {x A)s . (2.4) 

Here cr(x) and ((x) are the Weierstrass a- and (^-functions with periods equal to 2ir and 
2 A. Thus, equation ( |2.3| ) has the unique solution 



Now we turn to the momentum map equation 

C(x + A) - e- ih A-\x - A)C(x - A)A(x - A) = 27riifZ<5(x), (2.6) 

where Z = Res x=0 C(x + A). 

By using a generic gauge transformation we can diagonalize the field A. Then equation 



£>D takes the form of eq. (|2.3| ) 



C'(x + A) - e- lh D- l C\x - A)D = 2mK'Z'5(x), (2.7) 

where 

A(x) = T(x - A)DT~ 1 (x + A), C{x) = T{x - A)C'(x)T-\x - A) 
for some T and Z' = Kes x= AC (x) . We also have 

K' = T~ 1 (0)KT(0) = T _1 (0)e <g> e'T(O) = f®v\ <f,v>=l 
i.e. / = T~ 1 (0)e and e*T(0) = u*. According to ( 2.5) we find 



ij 



Taking the residue of C'(x) at the point x = A we arrive at the compatibility condition 

Z' = K'Z' = f <g> v*Z', <f,v>=l. 



The solution of this equation is Z' = f <g> g l , where g is an arbitrary vector. Now it is 
easy to find Z: 

Z = T(0)Z'T- 1 (0) = T(0)/ ® ^T-^O) = e <g> ^T^O) = e®b l 

Thus, we get 

C(x + A) - e~ ih A~\x - A)C(x - A)A(x - A) = 27rz(e ® e*)(e ® 6*)<J(ar), (2.8) 

where e (g) 6* is a residue of C(a;) at x — A. 

To summarize, eq. ( |2.6| ) has a solution for any field A and for any field C, having a 
residue at x = A of the form e (g) 6*. For a fixed field A and a vector 6 this solution is 
unique. Note that, in general, < b,e >^ 1. The form of the r.h.s. of ( |2.8| ) shows that the 
isotropy group of this equation is 

G isot = {T(x) C G(x) | T(0)e = Ae, A G C}. 

This group transforms a solution of ( |2.8|) into another one, so the reduced phase space is 
defined as 

all solutions of ( p.6| ) 



V, 



red 



Gisot 

Since the group G isot is large enough to diagonalize the field A, we can parametrize the 
reduced phase space by the section (D,L), where L is a solution of ( |2.6| ) with A = D. 
One can easily see that V re( i is finite dimensional and it's dimension is equal to 2N ', i.e. 
iV coordinates of D plus iV coordinates of the vector b. Due to eq. ( |2.5| ) the corresponding 
L-operator has the following form: 

Multiplying L(x) by the function ' 7 ^~^2'+h) e ^ ( x ~ A ) h i performing the gauge transforma- 
tion by means of the diagonal matrix e^~( x ~ A ^ q , and making the shift x — > x + A we 
obtain the L-operator of the elliptic Ruijsenaars-Schneider model: 

L n 3 (x) = e * e * q L(x + A)e *■ X! ? (2.10) 

<j{x + n.) 

Let us briefly discuss the Hamiltonian. It is well known that the simplest nontrivial 
Hamiltonian invariant with respect to the action of the current group is given by: 



H = I dxtrC(x), (2.11) 
where a is a constant. It is not difficult to show that on the reduced phase space 



2m N 

H ^=- mF , (2-12) 

i=l 



Nil - e~ ih 1 



that is up to a constant nothing but the simplest Hamiltonian of the elliptic Ruijsenaars- 
Schneider model. 
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3 The Poisson structure on the reduced space 



In this section we are going to prove that the Poisson structure on the reduced phase space 
does coincide with the Poisson structure of the elliptic Ruijsenaars-Schneider model. In 
other words we need Poisson brackets for the coordinates D-s and b-s. According to 
the general Dirac construction one should find a gauge invariant extension (we mean the 
invariance under the action of Gi SOt ) of functions on the reduced phase space V re d to a 
vicinity of V re d and then calculate the Dirac bracket. 

One can easily write down the gauge invariant extension for the matrices D and L(x) 
while the bracket for the coordinates -D, and b{ can be extracted from the bracket for D 
and L{x). This extension looks as follows: 

D -> D[A] = T~ x [A\{x - A)A{x)T[A\{x + A) (3.13) 

L{x) -> C[A,C](x) = T _1 [A](x — A)C{x)T[A\{x — A). (3.14) 

Some comments are in order. Eq. ( |3.13j ) is a solution of the factorization problem for 
A(x). Generally this solution is not unique but we fix the matrix T[A] by the boundary 
condition TL4](0)e = e that kills the ambiguity and makes ( p,13| ) to be correctly defined. 
It is obvious that on V re d- T[A] = 1 and C[A, C](x) = L(x). 

We start with the calculation of the Poisson bracket for C{x) and C(y). We omit the 
discussion of the contribution from the second class constraints to the Dirac bracket till 
the end of the section. By definition, one has 

{Ci,C 2 } Vred = ({T X ,T 2 }L X L 2 -L 2 {T X ,T 2 }L X -L X {T X: T 2 }L 2 

+ L X L 2 {T X , T 2 } + {d, C 2 } - {Tx, C 2 }Lx - {C x , T 2 }L 2 

+ L 2 {C x ,T 2 } + L x {T x ,C 2 })\ Vred (3.15) 

Here we took into account that T[A]|-p red = 1. 
Let us first calculate 



{Cij(x),T kl (y)} = Y^ J dz{C ij (x),A mn (z)} 



$T kl (y) 



SA mn [z) 

Performing the variation of the both sides of (|3.13|) , we get 

X(x) = t(x- A)D-Dt(x + A) +d, (3.16) 

where X{x) = SA(x), t(x) = 8T(x) and d = 5D. 
The general solution of ( |3.16| ) is 

1 r 1 

t(x) = Q - — ^2 J dz—w(x - z, q ij )X ij (z)E ij . (3.17) 

Here Q is some constant diagonal matrix and the function w(x, 0) should be understood 

as 

w(x, 0) = ]im(w(x, e) - - ) = ((x - A) - ^-(x - A) - ~. 

e^O 1 — e ie 7T Z 
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Note that these functions solve the equations 

^—.{w(x + A, qij ) - e~ iqij w(x - A, qij )) = S(x) - —Si,. 
2m 2tt 

The solution t(x) obeying the condition t(0)e = has the following form 

1 r 1 1 

K x ) = Y~ 5 / dz (jy w (- z , qij)X ij (z)E ii - —w{x - z, q ij )X ij (z)E ij ) (3.18) 

Performing the variation of eq.( p.!8| ) with respect to X mn {z) one gets 
ST kl (x) _ kl , . , 11 



V red = Q mn ( X + A > Z ) = n — : 7T~ ( W (- Z , Qkn)SklSkm ~ w(x ~ Z, q k l)S km Sl n ) 



5A mn (z)"- d ' 2mD,_ 

Thus on the reduced space we obtain 

-{C 1 (x),T 2 (y - A)}| red = Ki 2 (x,y)Li(x) - Li(x)ui 2 (x,y), 

7 

where 

Ki2(x, y) = tr 3 J dz(D 3 r™(x -z-2A)- r l *{x - z)D 3 )Q 23 (y, z), 
to 12 {x, y) = tr 3 J dz(D 3 r l *(x - z - 2 A) - r^{x - z)D 3 )Q 23 (y, z). 
We also get 

-{T x {x - A), C 2 (y)} = -Pk 12 ( V , x)PL 2 (y) + L 2 (y)Pu 12 (y, x)P. 

7 

By using the relation 

DjQtjix, z) - DiQ^x, z-2A) = S(x - z)S ik S jl - 5(z - A)S ik S kl = S$, 

we find 

Kij ki (x, y) = -r+ (x - y)ij M + ^ r + (x - A)y km 8 M 

m 

Uij u{x, y) = kij kl (x, y) + 2nD i Qf i (y } x). 

Recall that 

-{Ci(x),C 2 (y)} = -r^x-y^^CM-C^C^r^x-y) 

7 

+ C 2 (y)r+(x - y)d(x) + d{x)r_{x - y)C 2 (y), 

Substituting {C,T}, {T,C} and {C,C} brackets into ( |3.15| ) we can rewrite the {£,£} 
bracket in the following form: 

-{Ci(x),C 2 (y)}\ red = -L l (x)L 2 (y)k + (x,y) - k'(x,y)L 1 (x)L 2 (y) 

7 

+ L 1 (x)s~(x,y)L 2 (y) + L 2 (y)s + (x,y)L 1 (x), (3.19) 



where 

k~(x,y) = r-(x-y) + K 12 (x,y)-PK 12 (y,x)P-{T 1 (x-A),T 2 (y-A)}, 

k + (x,y) = r + (x - y) + uj 12 (x,y) - Puj 12 (y,x)P - {T x {x - A),T 2 (y - A)}, 

s~(x,y) = r-(x-y)+u; 12 (x,y)-PK, 12 (y,x)P-{T 1 (x-A),T 2 (y-A)}, 

s + (x,y) = r + (x-y) + K, 12 (x,y)-Pu> 12 (y,x)P-{T 1 (x-A),T 2 (y-A)}. 

It is easy to find P/c ± (x, y)P = —P5(x — y) — fc ± (y, x) and Ps ± (x, y)P = ±s T (y, x). We 
also have one more important identity 

k + (x,y) + k~(x,y) = s + (x,y) + s~(x,y). 

To complete the calculation we should find the bracket {Tij(x — A),T k i(y — A)} on 
the reduced space. The straightforward manipulations lead to a divergent result. By this 
reason we define this bracket as follows: 

{T l3 (x - A),T kl (y -A)} = \ Jim ({T y (x - A),T^(y - A)} + {T^x - A),T kl (y - A)}) 

where T^(x) is defined as a solution of the factorization problem with the boundary 
condition T(e)e = e. We have 

{T l3 {x - A),T^(y - A)} = J dztiQftJx, z)Q% e (y, z'){A mn {z), A sp (z')} 

= 7 / dzdz' (-r + (z - z') mn sp (D n Q% n (x, z) - D m QZ n (x, z - 2A))S% £ (y, z') 

-2nP mn sp 5(z -z' + 2A)D m Q% n (x, z)S% %, z')) 

One can prove the cancellation of the singularities as e — > 0. The result for the bracket 
{T, T} is 

-{Tij{x- A),T kl (y- A)} 
7 

= -r+(x - y)ij k i + J2 r +( X - kmhl + J2 r ~( A ~ V)im kAj 

m m 

1 11 

+-w(x -y + A, q ik )5 jk 5ii - -w(x, q ik )5 jk 5 k i + -w(y, q^S^Su 

% % % 

+-Sij5 ik 5u + -(C(qi k ) — ^^—qik)Sij5 k i(l — 5 ik ) — - y2(E ab — E ba ) ik 5ij5 k i 

^ 2 a<b 

Combining all the pieces together and taking into account the identity e~ ts w(x, s) = 
—w(—x, —s) we get the following expression for the coefficients: 

1 1 

k ij ki( x , V) = --C(Qik)SijSki{l - S ik ) - - (((x -y) + C(y - A) - C(x - A)) SijSikSu 

% % 

- - (w(x -y + A, q ik )5a5 jk + w(y, q^SuS^ - w(x, q ik )S jk 5 k i) (1 - 8 ik ) 

% 
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+ — — - ' Qik^ijSki + - y2(E a b — Eba)ik5ij5ki 
i tt 2 

a<b 

Kj ki( x > v) = \ (c( x -y)- ~^( x ~ y)j dijtiikSu 

+ - (w(x - y + A, q ik )5 jk 5u - C(fe)<%<^) (1 - 8 ik ) 

+ T ~ - Qik^ij^kl + 77 y2(E a b — E ba ) ik 5ij5kl 

I TT 2 ~i 

a<b 

$ij ki( x > y) = -\ [c{y - A) - ^-(y - A)^J s^Su 

- t (w(y, qki)SijS u + ((qi k )5ij8 kl ) (1 - 5 ik ) 

+ T ~ ^ Qik^ij^kl + - y2(E a b — Eba)ik5ij5 k l 

% TT 2 ~1 

a<b 

s tj ki( x > y) = \ (cO - A) - -^(x - A)^J SijSikSu 
+ - (w(x, q ik )S jk 5 k i - ((q ik )5 ij 5 kl ) (1 - 8 ik ) 

+ - qik^ijSkl + 77 y2(E a b — Eba)ik5ij5 k l 

% it 2 *-ri s 

a<b 

It is instructive to note that one can check by direct calculation that the term ^^-qi k 5ij8 k i 
in the expressions obtained for k-s and s-s does not contribute to the bracket {£, £}. 
Remind that (see eq.([2~^)): 

La{x) = -j=w{x,h)bi, (3.20) 

so to obtain the bracket {hi, bj} it is sufficient to examine the {La, Ljj} bracket only. The 
crucial point which can be checked by the direct calculation is that the bracket of Cu with 
the constraint ( |2.2|) vanishes on V re d in the limit e — > 0. Thus, there is no contribution 
from the Dirac term to the {La, Ljj} bracket. 

By substituting the expressions obtained above for k and s into eq. fl3.19Q , one gets for 

-{L u (x),Ljj(y)} = -Lji(x)L ij (y)w(x-y+A,q ij )--Li j (x)L j i(y)w(x-y+A,q j i). (3.21) 
7 i i 

It follows from this equation that 
*r, h \_ hh s ji) w (yi Sjj)w(x -y + A, qjj) - w(x, Sjj)w(y, Sjj)w(x - y + A, gfo) 

|Oj, Oj J — OiOj . . . . 

7 w(x,h)w(y,h) 

(3.22) 
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By using one of the known elliptic identities we get 



% -{b u b 3 } = Mj(2Cfe) - C(Qij + h)- ((% - h)). (3.23) 

To complete the examination of the Poisson structure on the reduced phase space one 
should find the bracket {£, D} and {D,D}. Performing the straightforward but rather 
tedious calculations following the same line as above, we find 

{D[A] 1 ,D[A] 2 }\red = 0, (3.24) 

-{£(x)x, D[A] 2 }\ red = - Lij{x) D 3 E ij ® %• (3-25) 

It is worthwhile to point out that there are no Dirac terms in these brackets because D[A] 
is invariant with respect to the action of the whole affine group GL(N). 

Now for the reader's convenience we list the Poisson brackets obtained in terms of the 
coordinates on V re d 

{qi,qj} = 
~{<li, b j} = bjdij 

-{Mi} = M,(2C(%)-C(% + ^)-C(%-^))- (3.26) 

7 



One can see that the dynamical system defined by (|3.26|) and ( 2.12 ) is nothing but the 
elliptic Ruijsenaars- Schneider model. 

In |7| the dynamical r-matrix for L Rui i was obtained by direct calculation with the 
help of the Poisson structure ( p.26|) . Comparing the r-matrix coefficients k and s with 
the ones in |?j we see that, in fact, they differ by the tensor ^J2 a <b(E a b — ^6o)ijfc%5*i- 
However, in our calculations of the bracket {£, £} we ignored the contribution from the 
Dirac term. We conjecture that just the Dirac term is responsible for cancelling this 
tensor. 



4 Equations of motion 

The equations of motion for the Hamiltonian ( |2.12j ) are given by 

D = {trL(x), D} = -"fL{x) diag D (4.27) 

and 

L(y) = {txL(x), L{y)} = [L(y), M(x, y)}, (4.28) 

where 

M(x, y) = -ji ^ (w(x, -q k i)L(x) k iE kk - w(x - y + A, -q M )L(x) k iE k i) . (4.29) 

kl 



2 It is interesting to note that this identity can be easily obtained from the x, y-independence condition 



for the r.h.s. of eq.( 3.22 ) 
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Here we use eq.( |3.19| ) and the explicit form of k and s. Since trL(x) is invariant function 



the contribution from the Dirac term vanishes. For the reader's convenience we note that 
by using the elliptic function identities || one can rewrite M = M(x + A, y + A) in the 
following form 

M = ll(x, h) ( C(x + Q-tt X ~ y) L(y + A) - (C(x + h)- C(a))E h)I (4.30) 

+ E^E(Cfe)-cfe-^))^-^E fo A ( 4 - 31 ) 

k i^-k k 

+ ^ aqu) -^ y + h h u {y + A)E H ), (4.32) 



where we have introduced l(x,h) = w(x + A,h). The first two terms in ( |4.3U ) are 



irrelavant, so M coincides with the standard M-matrix of the elliptic RS system. 

We show that the general solution of the equations of motion for the elliptic Ruijsenaars- 
Shneider model is given by 

D(t) = D[e- 2n ^ Lo{x)t D }, (4.33) 
where D = D[A] denotes the solution of the factoriztion problem ( |3.13|) : 



A(x) =T(x- A)D[A]T(x + A)- 1 (4.34) 

and D , and L (x) are the coordinates and the L-operator at t — respectively. 
To prove ( |4.33| ) we start with calculating the derivative D(t): 

where A t (x) = e~ 2n " ,L ° ( - x ^ t D . One can find the derivative by performing the varia- 

tion of eq.([0§: 

(T-\x - A)5T{x - A))D[A] - D[A]{T-\x + A)5T{x + A)) + 5D = X(x), (4.36) 

where the notation X(x) = T~ 1 {x — A)5A(x)T(x + A) was introduced. In contrast to 
( |3.16|) in eq. (|4.36|) we do not impose the constraint T — 1. 



Now we solve ( f4.36| ) for SD: 



dx 

5D = I —X(x) kk E kk . (4.37) 



Eq. fl4.36D also allows one to find the matrix 

f dz ( 1 1 \ 

T~ x {x)bT{x) = V / — —w(-z, q k i)X(z) M E kk - —w(x - z, q k i)X{z) k iE M 

2m \D k D k J 

(4.38) 

that will be used in the sequel. From Q4.37|) , (|4.38|) we find 



W[AIU 1 T£(z - A)T jk (z + A) (4.39) 
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and 



= ^^^T^-A^z + A) (4.40) 



1 w(x-z,q k i) x , 

~ T ki {z- A)T jl {z + A). 



2ni D k 

Substituting ( |4.40|) in eq. ( |4.35| ) and taking into account A t (x) = —27vyL (x)A t (x), we get 



D(t) kk = -7/ dzT^(z - A)L {z) im T mn {z - A)T^{z - A)A t (z) sj T jk (z + A) (4.41) 
that with the help of ( |4.34j ) reads as follows 



D(t) = -7/ dz(T-\z - A)L (z)T(z - A)) diag D(t). (4.42) 

The last formula implies the notation 

L t (x) = T-\x - A)(t)L (x)T(x - A)(t) (4.43) 

that provides the Lax representation ^L t (x) = [L t (x),M(x)] with M(x) = T~ 1 (x — 
A)T(x-A). 

Let us show that the Lax operator L t (x) coincides with the L-operator of the elliptic 
Ruijsensaars-Shneider model. To this end we calculate explicitly M(x). We have 



M kl (x) = J (t-\x-A) 



ST(x-A)\ , dA t {z Hj 



Substituting ( [4.40| ) and using the relation e ts w(x, s) = —w(—x, —s) we get 



M(x) = -ji dz^, (w{z, -q k i)L t (z) k iE kk - w(z - x + A, -q H )L t (z) k iE k ^ . (4.44) 

J kl 

Note that this expression literally coincides with ([4.29|) if we change L t for L. Since at 
t = the operators L and L are equal to L , they coincide for any t. 

5 Conclusion 

We have proved that the elliptic Ruijsenaars-Schneider model can be obtained by means 
of the reduction procedure. It is worthwhile to point out that we have used not the 
Hamiltonian but the Poissonian reduction technique. Our construction is specified by 
the choice of the trigonometric r-matrix for the Poisson structure on HD and by fixing 
the special value of the momentum map. By varying the r.h.s. of the momentum map 
equation one can derive some other systems. For instance, it is not difficult to specify the 
momentum map equation in a way that leads to the elliptic Calogero-Moser model. It 
clarifies the coincidence of the dynamical r-matrices for these two models pointed out in 
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We have considered the simplest example of HD for GL(N). It seems to be interesting 
to examine the Poissonian reductions of HD that correspond to some other choices of Lie 
groups or r-matrices. 
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